The imminent detection of gravitational waves from 
massive black-hole binaries with pulsar timing arrays 
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r] ! Recent observations of massive galaxies indicate that they double in mass and 
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quintuple in size between redshift z = 1 and the present, despite undergo- 
ing very little star formation, suggesting that galaxy mergers drive the evolu- 
tion. Since these galaxies will contain supermassive black holes, this suggests 
a larger black hole merger rate, and therefore a larger gravitational-wave sig- 
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^ \ nal, than previously expected. We calculate the merger-driven evolution of the 

mass function, and find that merger rates are 10-30 times higher and grav- 
; itational waves are 3-5 times stronger than previously estimated, so that the 

7— I ■ 

gravitational-wave signal may already be detectable with existing data from 
" pulsar timing arrays. We also provide an explanation for the disagreement 

£ : 

with past estimates that were based on dark matter halo simulations. 

Until recently, the massive galaxies living near the centers of clusters were thought to be 
"red and dead" during the interval z < 1 (equivalently, for D L < 6.7 Gpc or the last 7.8 billion 
years, assuming WMAP's seven-year cosmological parameters (7)), meaning they no longer 
underwent much star formation, and they remained essentially static. However, a growing body 
of evidence, both observational (2-4) and computational (5, 6), indicates that while this epoch 
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indeed saw little star formation in massive galaxies, these galaxies were far from static: to the 
contrary, massive galaxies are seen to have typically doubled their mass, and quintupled their 
size during this period (3). The most likely explanation for the observations is that massive 
galaxies merge with less massive satellite galaxies, and this is supported by simulations. Since 
these merging galaxies will generally also contain massive black holes, the larger-than-expected 
galaxy merger rate also implies an unexpectedly large black-hole binary merger rate. This has 
significant implications for pulsar timing arrays (PTAs) that are currently searching for these 
gravitational- wave sources. 

We therefore calculate the expected gravitational- wave signal from this entire population. 
Whereas the most recent estimates of this signal rely on dark matter halo simulations (e. g. (7)) 
to provide a population of merging binaries, we instead develop a simple observationally- 
normalized model for merger-driven galaxy evolution. We expect our resulting signal estimate 
to be significantly larger for the following reason: if we treat the satellite galaxy as a stellar 
bulge inside a halo, rather than a halo on its own, this has a surprisingly large impact on the 
mass of the satellite that survives being stripped by the host, and therefore the dynamical fric- 
tion timescale t di for the satellite to reach the center of the cluster. If we use an NFW profile (8) 
for the halo, it is simple to show that t di oc (1 + M*/M DM )~ 9 (9), where M* is the stellar mass, 
and M D m is the mass of the tidally limited subhalo containing the stellar bulge with an embed- 
ded black hole. Therefore, if after tidal stripping we find that M*/M DM « 1/2, for example, 
then the satellite will sink ~ 40 times faster than it would in a pure dark matter simulation, 
suggesting such simulations may seriously underestimate merger rates. 

We begin describing our model by noting that under the assumption of negligible star for- 
mation, the probability of a merger between black holes of mass M' and M" during the redshift 
interval z-z + dz is given by 

dJ ) dM"dz dM ' dM " dz = Pdz 0(M '' z) dM ' 0(M "' z) dM " ' (1) 



where is the black-hole mass function used in (10) and P is a normalization constant. We 
also note that the probability of finding a binary in a given frequency interval df is found by 
observing that the "chirp rate" / oc f 11 ^ 3 , so that ^ oc /~ 11//3 . We can therefore express 
the expected signal strength from the full population of binaries with respect to the probability 
distribution as 

pi /.lO n M pM 2 j4 

hl(f) = dz dM 2 dM 1 Nh 2 —-J-— , (2) 
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where N is the total number of sources, and h is the strain from an individual source. We im- 
plement Eq. © through Monte Carlo simulations, summing the power of each selected binary 
to arrive at the full stochastic signal power, h 2 c signal . The instrumental error from the radiometer 
is presently the limiting noise source for pulsar timing arrays, and places a limit on the amount 
of strain, h Ct no i SC , that is detectable in the data. The signal-to-noise ratio (SNR), p, is simply 
given by 
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where the strain signal is averaged over the orientation angles of the binaries and the sensitivity 
is averaged over the sky position angles. In Fig. [Q we show the signal strain for our present 
model and for the current standard result used in PTA analyses, as well as the noise level for the 
current generation of PTAs assuming a 5 year observation. We find that our mean SNR is ~ 4 
times the previous estimate; the triangles AACD and A ABE show the square root of the p 2 
integrand, so their areas are visual indicators of the relative SNRs of our signal and the previous 
estimate, respectively. Our mean signal estimate is ~ 5 x stronger than the previous estimate at 
/ = 1 year -1 ; indeed, our mean estimate of h Q > 5.8 x 1CT 15 , where h c = h Q (//l year -1 ) 2 ^ 3 
defines h , nearly equals the 95% null confidence limit of the European PTA (6 x 1CT 15 ) (11) 
and the NANOGrav PTA (7 x 1(T 15 ) (12). This indicates that some portion of the space of 
possible observational parameters is already ruled out within our model, and that we expect, in 



the statistical sense, that a positive detection will be made with PTAs in the very near future. 
If we assume that radiometer noise continues to dominate over red timing noise, then the noise 
floor of PTAs will improve with observation time T as T~ 13 / 6 (12). Therefore, by 2016, PTAs 
will be sensitive to our most pessimistic set of model parameters which yield h Q = 2 x 1CT 15 , 
so that we expect a detection by 2016 with 95% confidence based on our merger-driven model, 
with improved data quality and the possible combination of data from different PTAs only 
serving to accelerate the discovery. Since Advanced LIGO will not begin science operation 
until 2015, and will not reach its design sensitivity until 2018-2019, it is therefore likely that 
PTAs will make the first direct detection of gravitational waves, a landmark discovery which 
will usher in a new era of observational gravitational- wave astronomy. 

References 

1. E. Komatsu, et al.,Astrophys. J. Supp. Ser. 192, 18 (2011). 

2. A. R. Robaina, et al., Astrophys. J. 719, 844 (2010). 

3. P. G. van Dokkum, et al, Astrophys. J. 709, 1018 (2010). 

4. I. Trujillo, I. Ferreras, I. G. de La Rosa, Mon. Not. R. Astron. Soc. 415, 3903 (201 1). 

5. T. Naab, P. H. Johansson, J. P. Ostriker 699, L178 (2009). 

6. L. Oser, J. P. Ostriker, T. Naab, P. H. Johansson, A. Burkert, Astrophys. J. 725, 2312 (2010). 

7. V. Springel, et al., Nature 435, 629 (2005). 

8. J. F. Navarro, C. S. Frenk, S. D. M. White, Astrophys. J. 462, 563 (1996). 

9. S. T. McWilliams, J. P. Ostriker, F. Pretorius, Gravitational waves and stalled satellites from 
massive galaxy mergers at z < 1 (2012). In preparation. 

4 



10. Y.-T. Lin, J. P. Ostriker, C. J. Miller, Astrophys. J. 715, 1486 (2010). 

1 1. R. van Haasteren, et al.,Mon. Not. R. Astron. Soc. 414, 31 17 (201 1). 

12. P. B. Demorest, et al. (2012). arXiv: 1201.6641 [astro-ph.CO]. 

13. A. Sesana, A. Vecchio, C. N. Colacino, Mon. Not. R. Astron. Soc. 390, 192 (2008). 

14. A. Sesana, A. Vecchio, Class. Quantum Grav. 27, 084016 (2010). 

15. We wish to thank Alberto Sesana and Priya Natarajan for helpful discussions. Correspon- 
dence and requests for materials should be addressed to S.T.M. at stmcwill@princeton.edu 



5 




-8.2 -8 -7.8 -7.6 -7.4 -7.2 -7 

log(f [Hz]) 

Figure 1: Gravitational wave strain and strain sensitivity for a 5 year observation with PTAs. 
The red dashed line is the approximate strain sensitivity for current PTAs (11), and the green 
dashed line shows the previous estimate for the stochastic signal strength that is currently in 
standard use for PTA analyses (13, 14). The dark blue solid line corresponds to our mean 
estimate for the stochastic signal strength, with the blue shaded region bound by thin solid 
blue lines showing our 95% confidence interval for this estimate, based on the observational 
uncertainties of our model parameters. The light blue (AACD) and cyan (AABE) shaded 
regions show the area corresponding to the square root of the p 2 integrand, to be integrated over 
logarithmic frequency intervals as in Eq. ©, for our expected SNR of 8, and the SNR of 2 
expected from previous estimates (13, 14), respectively. 
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